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Abstract 

A patched-grid algorithm for the analysis of com- 
plex configurations with an implicit, upwind-biased 
Navier-Stokes solver is presented. Results from both 
a spatial-flux and a time-flux conservation approach to 
patching across zonal boundaries are presented. A gen- 
eralized coordinate transformation with a biquadratic 
geometric element is used at the zonal interface in order 
to treat highly stretched viscous grids and arbitrarily- 
shaped zonal boundaries. Applications are made to the 
F-18 forebody-strake configuration at subsonic, high- 
alpha conditions. Computed surface flow patterns com- 
pare well with ground-based and flight-test results; the 
large effect of Reynolds number on the forebody flow- 
field is shown. 


«i-e 

Nomenclature 
coefficients in surface definition 

e 

mean aerodynamic chord 

F, G, H 

flux vectors (convective and pressure terms) 

H v 

flux vector (viscous terms) 

Moo 

freestream Mach number 

Q 

time flux vector 

f 

position vector 

Rc 

Reynolds number based on c 

t 

time 


Cartesian coordinate system 

y + 

turbulent inner-law variable 

a 

angle of attack, deg 


body-fitted coordinates 
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Introduction 

There is an increasing effort in the development 
and application of Euler /Navier-Stokes computational 
algorithms for realistic aircraft configurations. Use of 
these prediction tools, in combination with ground- 
based and flight experiments, can significantly reduce 
the required time and cost of the design cycle. In 
the case of the recently proposed National AeroSpace 
Plane, a greater reliance must be placed on computa- 
tional methods since many of the high enthalpy and 
high Mach number conditions arc impossible to simu- 
late with present ground-based facilities. Timeliness is 
a critical issue in the design cycle and there is a pressing 
need to reduce the overall time required for an analy- 
sis of candidate configurations which might arise in a 
typical design cycle. 

The solvers in use today can be loosely classified 
into three types according to the underlying grid topol- 
ogy ; (1) Cartesian (2) unstructured and (3) multiblock 
structured. The Cartesian grid solvers lead to the sim- 
plest algorithms with the fewest operations and small- 
est memory requirements per grid point, but are least 
efficient and accurate in applications to curved sur- 
faces. With the latter two approaches, the use of body- 
conforming meshes leads to straightforward treatment 
and accurate resolution of general surfaces. The un- 
structured mesh solvers are the most general algorithms 
and their principal advantage is that grids over complex 
geometries can be generated in the shortest time period. 
However, these algorithms are also the most costly in 
terms of operations per time step and memory per grid 
point because of non-sequential memory access and the 
explicit connections required between the cells of the 
mesh. 

The multiblock structured meshes are more effi- 
cient in terms of operation count and memory than the 
unstructured meshes, since general connection informa- 
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tion is restricted to the edges of the blocks. However, 
more constraints are placed on the grid generation pro- 
cess in each block, although there is an active body of 
ongoing research to reduce or eliminate many of these 
constraints. The structured mesh approach leads to 
the most efficient algorithms for treating viscous flows, 
as the grids can be highly stretched in the direction 
normal to the developing shear layers. Of course, the 
algorithm of the future is likely to be a hybrid one, with 
overlapping structured meshes near the body linked to 
a Cartesian grid in the far field through an unstructured 
mesh, for instance. 

The present work is directed toward the develop- 
ment of improved techniques associated with the block 
structured approach, and in particular the patched grid 
approach. In this method, also referred to as a do- 
main decomposition method, the computational do- 
main is divided into a number of zones or blocks, each 
of which models the local geometric features and requi- 
site physics of the configuration. 1-8 The grids in each 
zone can be determined independently, with the net re- 
sult that the grid generation task is simpler and the 
solution is more efficient, since local clustering of the 
grid cells to resolve geometric and physical features in 
one region need not propagate to other blocks. The 
utility of the approach in a longitudinally-patched grid 
framework for an F-18 and an SR-71 configuration at 
supersonic speeds has been demonstrated previously. 3,3 
The method is extended in the present work so that 
the patching plane can be a general three-dimensional 
surface. The extension also overcomes a problem en- 
countered in the application to viscous flows, in which 
a highly stretched grid is present on either side of a 
zonal interface. 

Two approaches are described for maintaining the 
accuracy of the solution across the zonal interface. Ap- 
plications using the two patching algorithms are made 
to the F-18 aircraft, at subsonic high-alpha conditions, 
as part of an ongoing high angle-of-attack research pro- 
gram being conducted by NASA. Computed surface 
flow patterns from solutions to the Navier-Stokes equa- 
tions at a — 30 deg are compared to ground-based ex- 
periments conducted at NASA Langley and to in-flight 
experiments conducted at NASA Ames-Dryden using 
the High Alpha Research Vehicle (HARV). 9 The results 
are used to better understand the low Reynolds number 
ground-based tests, which are laminar and transitional, 
and their relationship to the flight results. A more com- 
prehensive comparison at a = 20 deg, including surface 
pressure comparisons, is given by Ghaffari et. al. 10 

Single Zone Algorithm 

The governing equations are the time-dependent 
Reynolds- averaged compressible Navier-Stokes equa- 


tions, cast in conservation law form and generalized 
coordinates as 


d§ dF 3G 

dt + d£ + dr} + 9$ 


( 1 ) 


expressing the conservation of mass, momentum, and 
energy. The superscript A denotes a quantity in gen- 
eralized coordinates. The thin-layer form of the equa- 
tions is used, where £ is the coordinate normal to the 
body surface. An ideal gas is assumed; the effect of 
turbulence is accounted for through the concepts of an 
eddy viscosity and eddy conductivity. The algebraic 
turbulence model of Baldwin and Lomax 11 is used to 
evaluate the turbulence quantities including the mod- 
ifications introduced by Degani-Schiff 13 to ensure the 
proper length scales are used in separated vortical flows. 

The equations, while written in generalized coordi- 
nates, are solved with a semi-discrete finite-volume al- 
gorithm, resulting in a consistent approximation to the 
conservation laws in integral form. 10,13,14 The convec- 
tive and pressure terms are differenced with the upwind 
flux-difference-splitting technique of Roe. 16 A MTJSCL 
(Monotone Upstream-centered System of Conservation 
Laws) approach of Van Leer 16 is used to determine 
state-variable interpolations at the cell interfaces. Each 
flux computation requires a single point on each side of 
the interface for first-order spatial differencing and two 
points on either side of the interface for the third-order 
upwind-biased spatial differencing used here. The shear 
stress and heat transfer terms are differenced centrally. 

The equations are advanced in time to the steady 
state in delta form, so that the steady state is inde- 
pendent of the time step. Two implicit algorithms are 
used :( 1) a spatially-factored diagonalized algorithm 13 
or (2) a hybrid algorithm, using streajnwise relaxation 
and approximate factorization in the cross-flow planes. 2 


Patching Algorithm 
General Concepts 

The two general approaches considered for patch- 
ing across zonal grids with coincident boundaries can 
be addressed based on the sketch in Fig. 1. Two- 
dimensional Cartesian grids in two zones are shown; the 
two zones have uniform spacing in the x-direetion and a 
grid mismatch in the y-direction. The indices (i, j) refer 
to the cell-center locations of zone 1 and (1, m ) to those 
of zone 2. The zonal interface across which the solution 
must be patched corresponds to the points defined as 
{Zt+l/ 2 ,i = 1 1 Jmax } and i m = 

Defining the spatial flux in the x-direction as F, 
Reu 1 has demonstrated that global conservation can be 
maintained by enforcing spatial-flux conservation along 
the interface, as 
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J F^ixi+trfdy = j F< 2) (sj_i /3 )dy. (2) 

The flux in zone 2, say, is constructed from the flux 
in zone 1 so that Eq. (2) is satisfied. For a scheme 
requiring two points on either side of the interface to 
compute the flux, the flux in tone 1 can be constructed 
from the data in zone 1 and an interpolation of zone 2 
data at a projection of zone 1 into zone 2. Referring 
to the conserved variables Q as the time flux of mass, 
momentum, and energy, much as F is referred to as 
the spatial flux of mass, momentum, and energy, the 
conservation of time flux in the region defined by the 
projected cells of zone 1 can be expressed as 

J J Q^dxdy= J J Q^dxdy. (3) 

It is assumed that the grid spacing normal to the inter- 
face boundary is nominally the same between the zones 
to be patched. If the spacing is not the same, the er- 
ror introduced is equivalent to that introduced into a 
single grid by a step discontinuity in the spacing. This 
approximation reduces the dimension of the interpola- 
tion implied by Eq. (3) by one, to a form very much liie 
Eq. (2). For equal spacing in the x-direction, Eq. (3) 
becomes 

J Q (1) (*<+i) d y= J <3 w (*i)dy. (4) 

On the other hand, the flux in zone I could have 
been constructed from zone 2 through Eq. (2), in which 
case a projection of zone 2 cells into zone 1 is re- 
quired. The time-flux constraint equation counterpart 
of Eq. (4) is then 

J Q^H x i~i) d y = J (5) 

Ideally, all of the Eqs. (2,4,5) would be satisfied. How- 
ever, any two of the three equations are sufficient to 
pass information between the zones. The approach 
which enforces Eq. (2) and either of Eqs. (4) or (5), 
as in the work of Rai 1,6 and Walters et. al, 2,3,8 is re- 
ferred to as the spatial-flux conservation approach. The 
approach which enforces Eqs. (4) and (5) is referred to 
as the time-flux conservation approach. This latter ap- 
proach satisfies Eq. (2) only to within truncation error 
but has been found, in practice, to maintain the con- 
servative, i.e. shock-capturing, properties of the single- 
zone scheme as well as the former approach. 7,8 The 
latter approach arises naturally in rezoning techniques 
for Lagrangian or adaptive mesh simulations as well as 
in transferring information between grids in a multigrid 


calculation. For higher order spatial differencing, the 
above equations need to be augmented with relations 
representing the additional information required at the 
interface. For a four-point stencil at the interface, such 
as used in the present work, a minimum of four inter- 
polations need to be done at an interface using either 
the spatial-flux or the time-flux appproach. 

Viewing Eq. (2) as a conservative interpolation of 
the time fluxes at a projection of one zone into an- 
other, it is clear that the approach does not rely on the 
boundaries being coincident; hence, overlapped grids, 
such as that used in the Chimera scheme of Benek 
et. al, 4 can be accomodated. Also, the time-flux ap- 
proach eliminates the departure from freestream con- 
ditions that can occur near curved interfaces using the 
spatial-flux approach, because of the differences in the 
discrete boundary definitions of the two zones. 8 

The conservative interpolations implied above can 
be evaluated discretely in a variety of ways. Defining 
the discrete flux in zone 1 to be interpolated at a par- 
ticular ae-position as cj x ^ from the discrete flux in zone 

2, Cm ' 1 , Rai 1 has used a piecewise-constant projection 
of C from one grid onto the other, as 


c l 1)= £ c m ) *r> (e) 

m 


N? = 


( 


0, 

0, 


if y m +i/ 2 < Vj-i/2 
if ym-l/2 > Vi +1/2 


u 


dy, 


otherwise. 


( 7 ) 

Nj* represents the relative area of overlap of cell m onto 
the cell j, 0 < N™ < 1, and the discrete flux balance 
is maintained, as 


E^w-E 0 *?**-. 

i m 

since = Ay m = y m +i /2 - y m -i/ 2 - The 

piecewise-constant approximation is, in general, only 
first order accurate. Interpolating onto a coarse grid 
from a fine grid, the coarse grid values represent an 
accurate integrated average of the fine grid solution. 
However, interpolating onto a fine grid from a coarse 
grid leads to the least accurate situation, since the dis- 
crete fine grid solution reflects the “stairstep* approx- 
imation of the coarse grid. This effect can be seen in 
the patched-grid inlet study of Ref. 7, where an oblique 
shock passes from a coarse grid into a much finer one. 
The virtue of the approximation is its simplicity and 
efficiency, especially in three-dimensional applications. 
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For second order accuracy, the distribution of flux 
in zone 2 can be taken as piecewise linear, 

C^{y) = <4 2) + {D v c&) m {y - *»), 


for < y < SWl/ 2 . (8) 


where the slope in each zone can be determined to 
maintain monotonicity between the zones, 16 as for ex- 
ample, 

/ 7(2) _ /M 2 ) 

(D„C^) m = Hr m ), 

!/m+ 1 — ym— 1 

(9) 

n 1 2r m 

<P( r m) - 2 

’ m » x 

(10) 

rm ~c^-c^ 

(11) 

The projection of the zone 1 grid onto the zone 
gives 

2 grid 

cj 1 > = J2 [ + (D y CW) m M? ], 

m 

(12) 


°» */ Vm+1/2 < 9J--1/3 

M r = | °> *7 s/m-i/2 > y,+i/2 

■t— / otherwise, 

(w) 

where the discrete flux balance is again maintained 
since M™Ayj — 0. This approximation repairs the 
inaccuracy of the piecewise constant integration when 
interpolating to finer grids from coarser grids. 

An alternate method within either approach is to 
use a nonconservative pointwise interpolation to evalu- 
ate the flux as from Eq. (8) evaluated at the cell centers 
of zone 2, 

cj l) = tf (2) (%) = c£> + (D v C™) m ( y, - y m ), 

for y m _i /2 < y i < y m +i/ 2 , (14) 

The discrete conservation constraint equations are then 
satisfied only to within truncation error. The lack of 
complete conservation is most apparent when interpo- 
lating to a coarser mesh from a finer mesh, where the 
flux interpolated at the center of the coarse grid zone 
may not be representative of an integrated variation 


of the finer grid flux over the coarser grid zone. This 
limitation can be removed by defining 

C 1‘’ - ^ £ C(,) M. («) 

p=l 

where C^(yj>) is defined pointwise from Eq. (8) eval- 
uated at 

2p — 1 

Vp = *fy-i/a + ~2 N~ Ay >’ P =1 >#> ( 16 ) 

representing the cell centers resulting from a division 
of the coarse cell zone into N smaller zones. This cor- 
rection need only be used when interpolating from fine 
cells onto coarse cells, so that N can be chosen auto- 
matically as, for example, 

N = max( i.JBilft-*- 1 /* ), ( 17 ) 

ym+l/2 — Vm- 1/2 

In the limit of IV — ► oo, the procedure satisfies the 
conservative property exactly. Exact satisfaction of the 
conservation equations is important for capturing weak 
solutions to the governing equations, such as shocks 
and slip surfaces. In smooth regions of the flowfield, 
Eq. (14) is sufficient. 

Surface Patching 

The counterpart in three dimensions of patching 
the zonal solutions along a line in two dimensions 
(Fig. 1) corresponds to patching along a surface. Given 
the three-dimensional coordinate transformation im- 
plied in Eq. (1) for each zone, 

(is) 

a patching surface can be defined without loss of gener- 
ality as a surface of constant £. The patching algorithm 
must match the solution between the zones given the 
set of surface points on each grid defining the patching 
surface, as illustrated in Fig. 2 for a planar interface. 

To interpolate across the interface, the transfor- 
mation defined by the discrete ordered set of points in 
zone 1, 

ffc) ■ * = l)*mo®5 3 = l>Jmaxi k" 1> lro<xsc}> 

is used to construct at the patch interface £ *= constant 
the generalised coordinates for the set of points in zone 
2 corresponding to the patch boundary, 

* rn = 1, rflfnax] n = 1, Sruu}* 
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The equation to conserve the spatial flux at the inter- 
face £= constant is 

J J fMdrtdi = J J P^dndi, (19) 

which is analagous to Eq. (2) for the two-dimensional 
case. Likewise, the counterpart to Eq. (3) is 


J J J Q^d(drid( = J J J 3 (1) d£di?df, (20) 


where the limits of integration span the cells adjacent 
to the interface. The spatial-flux conservation approach 
uses Eqs. (19-20) to patch across the boundary; the 
time-flux conservation approach uses Eq. (20) and an 
analagous equation for the projection of zone 1 into 
zone 2. For higher order differencing, Eq. (20) is aug- 
mented by a similar equation relating the conserved 
variables in the region spanned by the second set of 
cells in the £ direction adjacent to the interface. The 
assumption of equal spacing across the interface reduces 
Eq. (20) to a two-dimensional equation. 

Eqs. (19-20) express the conservation relations for 
the redistribution of flux on one side of the interface 
onto the cell faces of the grid opposite the interface. 
Assuming a piecewise constant variation of the flux, the 
interpolation reduces to determining the area of overlap 
from the cells (m, n) of one zone onto the cells 
(j, k) of the other, i.e. 




( 21 ) 


line segments of the two meshes, since all of the ar- 
eas of overlap are formed by the intersection of mesh 
lines from the two surfaces. This procedure has been 
used in the three-dimensional calculations of Kathong 
et. al,® Walters et. al, 3 and Ghaffari et. al. 10 The com- 
putational work for the area of overlap scales linearly 
with the number of mesh points. In a straightforward 
application, this procedure is more efficient than the 
clipping approach. 3 

The procedure used in most of the calculations re- 
ported here is to interpolate to the cell centers of one 
grid assuming a linear variation of the flux within the 
cells of the other grid, as 

c i,k — ^m,n + ~ >?m) 

+ (-D f C (2) )„,,„(& - £ n ), 


f or *1m— 1/2 ^ Vj — fn— 1/3 — Sfc — fn+l/2> 

analagously to Eq. (14) in the two dimensional case. 
For both approaches, the geometric information de- 
scribing the interpolation from one zone to another is 
calculated initially and then re-used at each subsequent 
iteration, so that the additional overhead due to the 
patching is minimal. 

In order to construct the generalized coordinates 
of one zone given the discrete transformation defined 
by another zone, a local geometric variation is assumed 
in each cell. A bilinear variation in each cell, 

fa + far) + fa$ + fat)t , (25) 


One approach to determining the area of overlap 
is a clipping algorithm, adapted from computer graph- 
ics, which has been used in the three-dimensional zonal 
calculations of Refs. 2,3, and 5. A more efficient pro- 
cedure, originally developed for Lagrangian hydrody- 
namic rezoning calculations by It amah aw, 17 relies on 
the divergence theorem applied to the position vector 
to calculate the area as a summation of line integrals 
over the bounding polygons. For a linear mapping of 
the surface, 

f= fa + far) + fat, (22) 

the area of a cell A p bounded by s directed line segments 
running from (*h,fi) to (*j 2 ,f 2 ) is 

3 ,- 1 m 

where ej is either +1 or —1 as the cell p lies to the 
left or right, respectively, of the line segment a. The 
areas of overlap can be calculated by sweeping over the 


leads to a mismatch in the discrete boundary definition 
between the two grids near a curved boundary, as illus- 
trated in Fig. 2. The boundary mismatch is not gen- 
erally an issue for meshes encountered in computations 
for the Euler equations, since the mismatch is a small 
fraction of the cure a of the cells adjacent to the bound- 
ary. For solutions to the Navier-stokes equations, in 
which the grids are highly clustered near the boundary 
to resolve the viscous layers at high Reynolds numbers, 
the situation depicted in Fig. 2(b) is not uncommon. In 
Fig. 3(a), the cell centers of one zone are shown in the 
generalized coordinate frame of the other zone (unit 
spacing in t) and f); the two cell centers nearest the 
boundary actually lie outside the discrete boundary of 
the other zone. For the zones shown, the grid spacing 
normal to the boundary is the same for both zones and, 
in the outer part of the grid, the cell centers lie midway 
between the grid lines, as expected. 

The surface definition can be substantially im- 
proved by using a degenerate biquadratic fit in the di- 
rection tangential to the boundary, i.e. 
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f= 3i+ 3317 + 23? + a 4 ij( 

+ ekt} 3 + 3ef/ 2 (. 

This form ensures (hat the boundaries of the cells are 
continuous between the zones, if the grid point data 
is augmented by additional edge points in the r) di- 
rection. For each cell edge along £ = constant, addi- 
tional points are determined by passing a least squares 
quadratic curve through the 2 grid points defining a 
cell edge and the nearest grid point on either side. For 
the zones in Fig. 2, the grid mismatch is reduced from 
0(1) to O(10~ 3 ), as shown in Fig. 3(b), wherein all the 
cell centers lie very close to midway between the grid 
points of the other. 

The generalized coordinates of each cell center to 
be interpolated is determined through a nonlinear itera- 
tive procedure. A nearest enclosing cell is guessed based 
on the distance from the grid point data to the cell cen- 
ter and then updated by solving for (t), f) from Eq. (26) 
through a Newton procedure, until the computed (rj, f) 
lies within the candidate cell boundary. The conserved 
variables are then interpolated using the linear fit of 
Eq. (24). 


Computational Results 
The grid topologies used are shown in Fig. 4. A 
2-zone grid, shown in Fig. 4(a), is described in detail 
by Ghaffari et. al. 10 The grids were generated in each 
block using an H-0 topology with transfinite interpola- 
tion. The two blocks are patched together at the strake 
apex. The upstream block contains 31 and 65 points in 
the circumferential and normal directions, respectively, 
and the downstream block contains 65 points in both 
directions. The 2-block grid contains a total of 185,000 
grid points and a spacing normal to the boundary of 
0.8 x 10~ 4 e. At the zonal interface, the grid lines are 
coincident only at the surface and in the longitudinal 
plane of symmetry and were generated with no require- 
ments for grid alignment in the interior of the grid. The 
disadvantage of the 2-block grid topology is that the 
circumferential grid clustering at the strake propagates 
forward into and ahead of the nose region, 

A three-block grid is shown in Fig. 4(b-c). Part 
of the upstream block in Fig. 4(a) is replaced with an 
0-0 topology. The additional block gives a much bet- 
ter resolution of the forebody geometry; the additional 
zonal interface is evident on the body surface and in the 
longitudinal plane of symmetry. A number of different 
grids were used in the first block. The finest, Fig. 4(b), 
has 65 and 73 points in the circumferential and normal 
directions, with a normal spacing of 0.1 x 10~ 4 C; the 
coarsest, Fig. 4(c), has 31 and 49 points in the circum- 
ferential and normal directions, with a normal spacing 
of 0.5 x IQ-*?. 


Laminar Flow Computations 

Symmetry plane pressure contours using the spatial- O' 
flux and time-flux conservation strategies for the 2- 
block grid are shown in Fig. 5 for the conditions con- 
sidered in detail by Ghaffari et. al 10 : = 0.6, R„ = 

0.8 x 10®, a — 20 deg. Both calculations interpolate 
across the zonal interface using Eqs. (21) and (23). The 
two solutions are nominally the same; any differences 
are on the order of the truncation error of the calcula- 
tions. 

Computed total pressure contours and surface par- 
ticle traces are shown in Fig. 6 at a high angle-of-attack, 
laminar flow condition: = 0.3, R c = 0.74 X 10 6 , a 

= 30 deg. The primary and secondary separation lines 
on the forebody are evident, as well as the secondary 
separation line on the strake. The vortices shed along 
these lines are very shallow and lie close to the body 
surface, especially in comparison to the primary vor- 
tex shed from the strake. The primary separation line 
on the forebody leads into a large region of reverse flow 
ahead of and predominantly under the strake apex. Un- 
der the strake along the body, a primary separation line 
is evident, leading to a vortex impinging on the strake 
lower surface. Along the fuselage above the strake, a 
separation line extending downstream from the apex is 
evident. The streamlines pass smoothly through the 
zonal interfaces. The overall flowfield is qualitatively 
similar to the laminar calculations given by Ghaffari ' v ~" / 
et. al 10 at ct = 20 deg. 

Oil flows from a test in the Basic Aerodynamics 
Research Tunnel (BART) 18 at NASA Langley Research 
Center for the same angle of attack but slightly lower 
Reynolds number, R c — 0.2 X 10®, are compared to 
the computed particle traces in Fig. 7. The experimen- 
tal oil flows demonstrate a striking resemblance to the 
computations. 

The effect of grid refinement is shown in Fig. 8 
for the 3-block grid at the experimental Reynolds num- 
ber. The interface interpolations were done with the 
time-flux approach, using Eq. (26) and the nonconser- 
vative interpolations of Eq. (24). With the finer grid, 
the forebody secondary separation line is straighter, ly- 
ing nearer to the leeward plane of symmetry, in better 
agreement with the experiment. The flow over the two 
downstream blocks is largely the same. Some apparent 
differences are due to the fact that the seeding positions 
for the particle traces are not identical for Fig. 8(a) and 
8(b). As can be seen by comparing Figs. 7 and 8, un- 
der the laminar assumption there is very little effect of 
Reynolds number over the range considered. In Fig. 8, 
the particle traces originating in a block are confined to 
remain in that block; the surface flow, however, passes 
smoothly through the two zonal interfaces. The first ' 
zonal interface has a variation in both the circumfer- 
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ential and normal spacing of the grid on either side, 
and the improved interpolation obtained with the bi- 
quadratic surface fit to the geometry is critical in ob- 
taining the smooth results shown. 

Turbulent Flow Computations 

Comparisons with wind-tunnel results at higher 
Reynolds numbers and with in-flight results show mark- 
edly different surface oil flow patterns. 9 The flight-test 
results correspond to a Reynolds number of 10-15 mil- 
lion based on S and show primary and secondary sepa- 
ration lines on the forebody much closer to the leeward 
plane of symmetry. Computed particle traces at flight 
Reynolds number for the 3-block grid of Fig. 4(b) are 
compared with the flight-test results 10 in Fig. 9. The 
primary and secondary separation lines on the fore- 
body are well predicted. The extent of separation in 
the strake apex region is considerably reduced and the 
secondary separation line on the strake is moved out- 
ward from the laminar condition, as expected. The 
computed secondary separation line on the strake doeB 
not extend into the forward region of the strake; a re- 
flex in the particle traces is indicated in that region, 
which generally leads with grid refinement to a clearly 
defined separation line. In this context, the grid spacing 
in block 1 corresponds to y + « 2 and is considerably 
finer in the normal direction than in either of the two 
downstream blocks. 

Concluding Remarks 

A patched grid algorithm for complex configure 
tions has been described. The work is an extension to 
the longitudinally-patched approach of Refs. 2-3, allow- 
ing for the analysis of grids which are highly stretched 
in the normal direction to resolve viscous flows, and 
for arbitrarily-shaped patch surfaces. Two algorithms, 
a spatial-flux and a time-flux conservation approach, 
have been used across zonal interfaces with few dif- 
ferences noted. The latter approach is somewhat more 
flexible and can be easily extended to handle more com- 
plex situations, such as overlapped and embedded grids. 
A longer term objective of the present domain decom- 
position approach is for the only constraint on the grid 
topology to be that the grids span the entire physical 
domain. The computational algorithm should be gen- 
eral enough to automatically determine the necessary 
connections between the domains to ensure a globally 
second-order accurate solution. 

Applications using the algorithm have been made 
to the F-18 forebody-strake configuration at a = 30 
deg over a range of Reynolds number. The laminar flow 
calculations agree well with the wind-tunnel results at 
Re — . 2 x 10®. The turbulent flow calculations are sub- 
stantially different in the separation line positions and 


agree well with flight-test results at R e = 10x10®. The 
wind-tunnel results at R e « 1 X 10® exhibit a separation 
pattern different from either calculation; 9 the difference 
is associated with the transition from laminar to tur- 
bulent flow. The laminar separation bubble described 
in Ref. 9 from flow visualizations occurs in the region 
of the primary separation line predicted by the present 
algorithm, assuming a laminar viscous model. 
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Zone 1 Zone 2 



Fig. 1.- Two-dimensional zonal interface. 
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(b) Nearfield view. 


Fig 2.- Coincident zonal interface defined by zone 1 grid (solid lines) and zone 2 grid (dotted lines) for partial 
ellipse section. 
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(a) Bilinear basis. 


(b) Biquadratic basis. 


Fig 3.- Generalised coordinates of zone 2 cell-center locations (dotted lines) in transformation defined by zone 1 
grid (solid lines) for partial ellipse section. 
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(a) Spatial- Bex conservation 
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(b) Time-flux conservation. 


Fig 5.- Centerline pressure contours. Moo — 0.6, /?,. — 0-8 x t0 r ’ f a = 20 cleg. 
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Fig 6,- Told pressure contours and surface particle traces for forebody-slrake using 3-bioclc grid. W* ---* 0J, R 
0.74 x 10*1 a ™ 30 deg. 
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(a) Side view. 

Fig 7.- Wind-tunnel (top) and computed (bottom) surface flow patterns on fore body-st rake at low-speed, laminar 
flow conditions, a — 30 deg. 
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(b) Three-quarter side view. 
Fig Concluded. 
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(a) 3-block grid {282,000 grid points). 

Fig 8." Effect of grid refinement on .forefeody-strake configuration. Side view, Af - 0.3, It,-. = 0.2 x i0 e a - 
30 deg. 
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(a) Side view. 


Fig 9.- Flight* test (top} and computed (bottom) surface tows for F-18 oafigmUtm. 
10®, a— 30 deg. 


ss 0.3, R e « Hi * 
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